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PACS numbers: 

In our entangled many-body quantum simulations we 
obtained three lines of evidence that mean-field-like dark 
solitons decay when taken as an initial condition: filling- 
in of the pair correlation function g^'^'; inelasticity of 
soliton-soliton collisions; and comparison of the exact 
quantum depletion of the macroscopic mode associated 
with the mean-field versus what one obtains with a 
weakly interacting Bogoliubov description. The latter 
is covered in a great deal more detail in [l|. Dziarmaga 
et al. have showed that fiUing-in of g^^^ does not nec- 
essarily correspond to filling-in or decay of a soliton in 
a single quantum measurement; that is, the soliton may 
delocalize, so that it fills in on average, but retain its 
soliton-like character in any single experiment, in keep- 
ing with well-established previous studies Q] ■ Martin and 
Ruostekoski have verified this result with an independent 
method 0,13. 

Dziarmaga et al.^s simulation of a single-shot measure- 
ment treats the near-mean-ficld regime in which the fill- 
ing factor for each site is 5000/31 ~ 161. They also 
work with an assumed mean- field soliton profile of a tanh 
function centered at different positions to build the de- 
localized soliton state. Our work uses a filling factor of 
1/2 to 2, and is not in the mean-field regime (although 
in the ground-state phase diagram our choices of hop- 
ping, interaction, and filling do correspond to the super- 
fluid region, not the Mott- insulating one) . We agree with 
Dziarmaga et al. that filling in of g^"^^ does not necessar- 
ily correspond to soliton decay in a single measurement. 
However, the common difficulty of connecting results of 
near-mean-field simulations to entangled quantum many- 
body simulations remains; large filling factors have a dif- 
ferent phenomenology than small filling factors, and one 
cannot use simulations from one regime to make a state- 
ment about simulations in the other. Dziarmaga et aZ.'s 
results are in a different regime from ours, and moreover 
have not dealt with our other two pieces of evidence for 
dark soliton decay. 

In a separate work |5|-l7|, one of us has studied the 
quantum analog of dark solitons over all interaction 
regimes, from a mean field Bosc-Einstcin condensate to a 
Tonks-Girardcau gas. We showed that it is in fact yrast 
states that properly play the role of phase slip, not in 


general mean-field-like dark solitons. These yrast states 
are precisely Lieb's type-II excitations, which were iden- 
tified with dark solitons in the weakly-interacting limit 
only. It is an intriguing open question as to whether or 
not anything like a dark soliton exists in the medium to 
strongly interacting regime, i.e., a well-defined, robust 
density notch which collides clastically with other soli- 
tons. In either case, the mean-field dark soliton is very 
far from the quantum dark soliton, as we have shown 
in [7|, and thus we believe that mean- field dark solitons, 
e.g., Dziarmaga et aL's superposition of tanh conden- 
sates, decay when injected into the quantum theory out 
of the near-mean-field regime (we note that the initial 
conditions in our Letter are close to symmetry-broken 
versions of Dziarmaga et al.'s tanh function). The spe- 
cific parameters for which mean-field theory breaks down 
for this problem were obtained via the Bethc ansatz, and 
the connection to the delocalization concept is discussed 
in more detail in [8|. 
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